Nonmagnetic impurity effects in MgB2 
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We study nonmagnetic impurity effects in MgB2 using tfie quasiclassical equations of superconduc- 
tivity for a weak-coupling two-band model. Parameters in the model are fixed so as to reproduce 
experiments on MgB2 as closely as possible. The quasiparticle density of states and the specific 
heat are calculated for various values of the interband impurity scattering. The density of states 
changes gradually from a two-gap structure into the conventional single-gap structure as the in- 
terband scattering increases. It is found that the excitation threshold is not a monotonic function 
of the interband scattering. Calculated results for the specific heat are in good agreements with 
experiments on samples after irradiation. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Superconductivity in MgB2 has attracted much at- 
tention since its discovery by Nagamatsu et alX Be- 
sides its high transition temperature = AQK for a 
phonon-mediated pairing mechanismyS it has a striking 
novel feature that approximately two different energy 
gaps open on different pieces of the Fermi surface. Ac- 
companied with theoretical predictions (.MiM, this multi- 
gap structure has been established within a short pe- 
riod by specific heat experimentSjiiSi^iiSiiiiiSii^ point- 
contact spectroscopyfi^ii^^ and scanning tunneling 
microscope^iiSiiS The larger gap opens on cylindrical 
Fermi surfaces of the cr-band with 2Ao-(0)/fcsTc — 3.6 — 
4.5, whereas the smaller one has the size 2IS.t,{Q) /ksTc ~ 
1.1 — 1.5 on three dimensional Fermi surfaces of the tt- 
band. 

Impurity effects of this material are very interesting 
to study. According to Anderson's theorem for classic 
s-wave superconductorSfSSi^ nonmagnetic impurities do 
not affect superconducting properties in zero magnetic 
field. However, it was shown later by Markowitz and 
KadanofS^ that is actually reduced in the presence of 
gap anisotropy and impurity scattering. An application 
to a two-band model is due to Golubov and Mazin.^'^ In- 
deed, they predicted a rather drastic decrease of Tc due to 
the interband impurity scattering. They also found that, 
as the interband scattering increases, the density of states 
changes from the two-gap structure inherent to the two- 
band model to the conventional single-gap structure. 
This reduction of Tc has been confirmed recently by a 
couple of experiments. Wang et ali^ measured the spe- 
cific heat of polycrystalline MgB2 after irradiation. They 
found both suppression of T^, and a tendency towards a 
single-gap structure as the scattering is increased by ir- 
radiation. Lee et alM clarified the possibility of com- 
plete suppression of superconductivity by replacing B in 
MgB2 by C. However, few quantitative calculations have 
been performed on the impurity effects based on a realis- 
tic model for MgB2.^^ For example, the specific heat has 
been calculated by Choi et alA and Golubov et alm^ based 
on the clean-limit Eliashberg equation to obtain excellent 



quantitative agreements. However, no detail study has 
been performed for the specific heat from clean to dirty 
limits based on a microscopic model. 

With these observation, we investigate nonmagnetic 
impurity effects in MgB2. We thereby clarify impurity- 
concentration dependence of the quasiparticle density of 
states and the specific heat, choosing the parameters in 
the model suitable for MgB2 . Section II gives the formu- 
lation. Section HI presents calculated results. Section IV 
summarizes the paper. We put fi = A;B = l throughout. 



II. FORMULATION 
A. Quasiclassical equations 

We start from the Eilenberger equations^ for the Suhl- 
Matthias- Walker modelSl with impurities, which form 
one of the most convenient frameworks to study impurity 
effects of the two-band model. We here adopt the formu- 
lation on the real energy axis instead of using Matsubara 
frequencies, which has an advantage that the free-energy 
functional can be defined unambiguously. 

The Eilenberger equation on the real energy axis is 
given for the uniform case by 

Here /^=/^(£, kp) and g^=g^ {e, kp) are retarded qua- 
siclassical Green's functions specified by the band in- 
dex a (= cr, tt) and the Fermi wave vector kp. They 
are connected by = (1 - f^f^'^f'^ with f^'^{tM) 
e, — kp). The symbol (• • • ) denotes the average 
over the Fermi surface for the relevant band with (1) = 1. 
The quantity r^^ is the relaxation time for the nonmag- 
netic impurity scattering from a- to /3-band; they satisfy 
7^ ~ ' ^'^^'^^ (0) normal-state density 

of states at the Fermi energy for the a-band. Finally, Aq, 
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is determined self-consistently by 



(2) 



where Sc is the cutoff energy, ip{e) =ta.nh{e /2T) , and Xa/s 
is dimensionless coupling constant with Aq^ = ^^Xpa- 
The functional for the free-energy difference corre- 
sponding to Eqs. (0 and (jJl are given by 
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Equation (PJ is a direct extension of Eilenberger's free- 
energy functional^^ to the two-band model. Indeed, vari- 
ations of Fs — Fn with respect to f^^ and A* lead to 
Eqs. (P) and lj2Jl, respectively. 

The entropy is obtained from this free-energy func- 
tional by Ss — Sn — d{Fs — Fn)/dT. Noting the station- 
arity of F^ — F„ with respect to /^^ and A* , we only have 
to differentiate with respect to the explicit temperature 
dependence in 0(e) = tanh(e/2T). We thereby obtain an 
explicit analytic expression for the entropy as 
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In deriving this expression, we have used Eq. ^ . Finally, 
the specific heat is calculated by numerically differenti- 
ating Eq. (|SJ) as 



Cs — C^n "I" T 



djSs - Sn) 

dT 
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Equations ^ and © form a convenient and efficient 
starting point to calculate the specific heat for various 
impurity concentrations. 

There is a disadvantage in the coupled self-consistency 
equations 1^ and ^ that they may not be very stable 
numerically. However, it can be removed when Sc is much 
larger than both A^ and 1/tq,^. This condition is well 
satisfied in MgB2 where £c corresponds to the Debye en- 
ergy ~ lOOOif iS Then using the asymptotic property 
^ Aa/ {—ie) as |£|— >cx3, Eq. (O is transformed as 
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where 7 = 0.577 is the Euler constant, £„ = {2n+l)TrT 
is the Matsubara frequency, and /„(£„) = f^{ien)- The 
limit Ec — > 00 has been taken safely in the first integral to 
transform the integration into the summation over Mat- 
subara frequencies. Equation Q tells us that A^^ has no 
angular dependence within the band. It then follows that 
fa neither has any angular dependence, so that (/„) — fa 
and {ga) — ga- Hence Eq. Q with e — > ie„ is simplified 
into 



£n + 
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Thus, the intraband scattering does not affect supercon- 
ductivity in zero magnetic field at all, in agreement with 
Anderson's theorem. Equation © could be presented 
at the beginning in place of Eq. However, Eqs. 

(Q)-® have an advantage that they could easily be ex- 
tended to nonuniform systems by simply adding gradient 
terms 

Equations JT)) and Q enable us to obtain A^ from 
calculations on Matsubara frequencies, which are numer- 
ically more stable than Eqs. and Once Aq, is 
fixed in this way, the density of states is calculated by 
solving Eq. Q as 

iV(e) = ^7V„(0)Reg:^(e). (9) 



B. Transition temperature 

We now derive the Tc equation valid at all impurity 
concentrations. To this end, we expand and ga up to 
first order in A^ as fa = fa"^ and = 1. Substituting 
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them into Eq. (|SJ, we obtain fa as 
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where hafs is defined by 
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Substituting Eq. (|10() into Eq. 0, we obtain the condi- 
tion for a nontrivial solution as 



det[l - i7] , 
where the matrix H is defined by 
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By solving Eq. Tc is obtained for an arbitrary t. 

Notice that all the summations in Eq. (|14|l can be ex- 
pressed in terms of the digamma function ipi^)- 

When TcT <C 1, Eq. l(T^ can be solved exphcitly by 
using the asymptotic expression ^(x) ^\nx {x~>-oo) as 



2e'' 

Tc = £cexp 

TT 



{Xapna-Xaanp)MecT) + l 



(15) 

with Aq = X^/jAq/j. This expression is useful to see 
whether is suppressed completely or not as the in- 
terband scattering increases. 



C. Density of states in the dirty limit 

We now derive an analytic expression for the density 
of states in the dirty limit of TcT <C 1. In this case with 
£<rc, we can neglect the first terms on both sides of Eq. 
(III). It then follows that /ffS /5^> or equivalently, 
/« ^ ff = and = = g"^. The quantities 



and g^ are obtained easily as 
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where e+ = e-|-iO+, and A is defined by 
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TABLE I: Coupling constants 
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Case A 


2.45 X 10"' 5.95 x 10~^ 


4.28 X 10"^ 


1.10 X 10"' 


Case B 


2.25 X 10"' 2.35 x 10"' 


1.69 X 10"' 


-1.22 X 10"' 



It hence follows that the density of states oc Re g^ has 
the conventional BCS structure, diverging at e = itA. 
The corresponding pair potential is determined by 



A„ = 1^ r &</>(e)/^(e) 



(18) 



as already pointed out by Golubov and Mazini2^ Thus, 
the ratio of the two pair potentials is given by Ao-ZA^r = 

Ao-/A,r- 



D. Numerical procedures 

There are five parameters in the model: 

Xaa, Xa^, K^, £c, A^^ (0)/A^^ (0) . (19) 

They are fixed so as to reproduce experiments on MgB2 
as closely as possible. To be more specific, the ratio 
N„{Q)/Nt^{Q) is set equal to 0.72 following an electronic- 
structure calculation. ^2 As for the coupling constants, we 
use Eq. (0) in the clean limit at T = 0, Tc which yields 



l = A..ln(g)+A..^ 
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/ 2e'''ec\ _ Ao-o-+ A7r7r + -\/(Acro-— A,r^)2 +4Acr^A^o- 
V ""^c / 2(ACTCTA,r^ — A(j^A7r(j) 

(20c) 

respectively. These equations are used to eliminate Xap 
in favor of Acro/A^o, Tc /A^q, and ec/Ao-o, where Ac,o,7rO 
denotes Ao-,7r(r = 0) in the clean limit. 

We here consider the following two cases i^SiiLi^ 

Case A: -^=0.50; ^ = 3.00; -^=20.0, 

ActO AttO Ao-O 

(21a) 



Case B : 



Tc_ 

A, 



0.48; ^ = 2.95; 



10.0. 



(21b) 

These values are chosen so as to reproduce tempera- 
ture dependence of the observed energy gaps on clean 
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FIG. 1: Transition temperature as a function of log j^Q(Ao-o''"CT7r) 



sainplesi^*iS4iSii& as closely as possible. The correspond- 
ing coupling constants are listed in Table I. It has been 
found that the whole results are rather insensitive to £c, 
as may be expected. 

Numerical calculations have been performed as follows. 
First, fa and ga with = 1 are expressed conve- 
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FIG. 2: The pair potentials as a function of T /Tao for three 
different impurity concentrations, (a) Case A; (b) Case B. 
Here Ta-r: is given in units of A~o 




FIG. 3: Density of states at T = O.OlTc for three different 
TcTTv- (a) Case A; (b) Case B. Here To-ir is given in units of 



nicntly in terms of a single function 
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Substituting Eq. and a trial (ActjAtt) into it, Eq. 

^ is transformed into a set of nonlinear equations for 
flc, which may be solved by using one of the standard 
numerical procedures.-'^ The obtained is then substi- 
tuted into Eq. ^ to find a new (Ao-jA^). This proce- 
dure is repeated until the convergence is reached. The 
pair potentials thereby obtained are then used in Eq. 
to calculate and on the real energy axis. Finally, 
those and g^ are substituted into Eqs. lO, JHl, and 
© to calculate the specific heat and the density of states. 



III. RESULTS 

Figure n plots Tc as a function of the interband scat- 
tering specified by logjo(^(TO "'"cttt) for the two cases of Eq. 
(|21|l . Here Tc is normalized by Tcq without the interband 
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FIG. 4: Specific heat as a function of T/Te for r^vr = 10000.0, 
1.0, and 0.1. (a) Case A; (b) Case B. Here Ta-n is given in 
units of A~n 




].2 0.4 0.6 0.8 1.0 1.2 
(a)° 




3.2 0.4 0.6 0.8 1.0 1.2 

T/T 



FIG. 5: Comparisons of theoretical curves with specific-heat 
measurements before irradiation (top), after the first irradi- 
ation (middle), and after the second irradiation (bottom) of 
Ref.Iia. (a) Case A with A<^oT-,T,r = 50.0, 7.0, 4.0 from top to 
bottom, (b) Case B Aao-r<T,r = 10000.0, 30.0, 1.0 from top to 
bottom. 



scattering. We observe that Tc drops steeply around 
T(j7r — A~Q^ in both cases. However, whereas Tc in Case A 
remains finite for Taj- ^0, in Case B decreases to zero 
at a finite To-Tr- This difference can be realized from Eq. 
(|15|l where the term in the square bracket passes through 
negative infinity for T^^O. It is interesting to see exper- 
imentally whether Tc of MgB2 is suppressed completely 
or not by increasing the interband scattering, although 
this may not be easyiS^ 

Figure|2]shows temperature dependence of the pair po- 
tentials for (a) Case A and (b) Case B. In both cases, the 
pair potentials decrease as Ta-n becomes shorter. How- 
ever, the two pair potentials do not approach to a single 
value even for Tc-k — * 0. Notice also that the pair po- 
tentials are not directly connected with any observable 
quantities except in the clean limit. 

Figure|21shows the density of states at T jT^ = 0.01 for 
(a) Case A and (b) Case B. In the clean limit A^oTaTr = 
10000.0, we clearly observe a couple of divergences at £ = 
Act, A^. As the interband scattering becomes larger, the 
divergences are smeared to finite peaks, which eventually 
merge into a single peak in the dirty limit at e = A given 



by Eq. (|17|l . Notice that the excitation threshold is not a 
monotonic function of Tcttt] for Case A, for example, the 
threshold is seen to increase as To-tt becomes smaller. 

Figure 01 plots the specific heat as a function of T/Tc 
in (a) Case A and (b) Case B for A^qt^t: = 10000.0, 
1.0, 0.1. A shoulder is clearly seen around T/Tc = 0.2 
for AcroTo-TT — 10000.0, corresponding to the divergence 
at e = Ajr in FigO It gradually disappears as To-^ be- 
comes shorter, however, and we finally have a conven- 
tional single-exponential behavior in the dirty limit. 

Finally, Fig. [S] compares the present theory with the 
specific-heat experiment after irradiation performed by 
Wang et alM^ The data points correspond to measure- 
ments before irradiation (top), after the first irradia- 
tion (middle), and after the second irradiation (bottom), 
where the former two are shifted upwards by 1.5C„ and 
0.75C„, respectively. The curves of Fig. 5(a) are obtained 
by using the parameters of Case A with A 50.0, 
7.0, 4.0 from top to bottom, respectively. On the other 
hand, those of Fig. 5(b) correspond to Case B with 
ActoTcttt = 10000.0, 30.0, 1.0. The agreements are good 
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for both cases, especially for Case B. From these com- 
parisons, we realize that the interband scattering af- 
ter the second irradiation is still not very strong with 
Ao-oT'tTir ^ 1-0. The fact may imply the difficulty of in- 
troducing the interband impurity scattering in MgB2.^^ 
This tendency is also seen in the carbon-substituted sys- 
tem Mg(Bi_a;Cx)2 where the two-gap structure is still 
observed clearly for x^O.lMi Thus, the large reduction 
of Tc observed in Mg(Bi_2;Ca;)2 by Lee et alm^ should 
be attributed not only to the interband scattering alone 
but also to a change of the pairing interaction due to the 
electronic structure. 



IV. CONCLUSION 

We have studied nonmagnetic impurity effects for 
MgB2 based on the quasiclassical equations of supercon- 
ductivity for the Suhl-Matthias- Walker model. The pa- 
rameters in the model are fixed so as to reproduce exper- 



imental values for A^o/^o-o and Ag-o /Tc. The interband 
impurity scattering tends to reduce the gap anisotropy. 
We have clarified how the density of states changes from 
the two-gap structure in the clean limit to the single-gap 
structure in the dirty limit with strong interband scatter- 
ing. Especially, there may be cases where the excitation 
threshold increases as the scattering becomes stronger. 
Calculated curves for the specific heat agree well with 
measurements before and after irradiation. This compar- 
ison has also enabled us to estimate the relaxation time 
TcTT for the interband scattering. It satisfies To-tt^I/Acto 
even after the second irradiation, implying the difficulty 
of introducing the interband scattering in this system. 
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